Preparation of Dicke States in an Ion Chain 
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We have investigated theoretically and experimentally a method for preparing Dicke states in 
trapped atomic ions. We consider a linear chain of A'' ion qubits that is prepared in a particular 
Fock state of motion, |m). The m phonons are removed by applying a laser pulse globally to the 
A'^ qubits, and converting the motional excitation to m flipped spins. The global nature of this 
pulse ensures that the m flipped spins are shared by all the target ions in a state that is a close 
approximation to the Dicke state jDjj"'). We calculate numerically the fidelity limits of the protocol 
and find small deviations from the ideal state for m = 1 and m = 2. We have demonstrated 
the basic features of this protocol by preparing the state |-D2^') in two ^^Mg"*" target ions trapped 
simultaneously with an ^^Al^ ancillary ion. 
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Entanglement is a fundamentally nonclassical feature 
of quantum mechanics and has been recognized as an im- 
portant resource in quantum information science [1] . The 
coherent manipulation of entangled quantum states can 
be useful for important tasks such as computation, cryp- 
tography, and metrology. In addition, the insights gained 
from research with entangled states may contribute to a 
better understanding of the foundations of quantum the- 
ory and shed light on physics at the boundary between 
the quantum and classical worlds. 

Various entangled states have been experimentally 
demonstrated in a wide variety of physical systems in- 
cluding photons, condensed matter systems, atoms in op- 
tical lattices, and trapped ions [2]. In the case of trapped 
ions, Bell states have been generated with high fidelities 
[3, 4]. For larger numbers of ion qubits, GHZ states [5, 6] 
and W states [7] have been realized. As the number of 
physical parts in a quantum system grows, we can find 
more classes of entangled states [8] and potentially more 
interesting applications, but the state preparation and 
verification becomes more challenging. Here we explore 
numerically an approach to preparing multiple trapped 
ions in a general class of entangled states called Dicke 
states. We then demonstrate the basic protocol on a 
mixed-species ion chain composed of two Mg+ ions cou- 
pled to a single A1+ ion. 

The Dicke state I-Dn"') is the equal superposition of all 
basis states of N qubits having exactly m excitations [9] . 
If we denote the qubit states as |J,) and |t), we can write 
an arbitrary Dicke state in the following form 
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The sum is over all (^) permutations (produced by the 
permutation operator Pk) with m qubits in the state jf). 
The W states are a special case of |D5v"0 with m = 1. 
Entanglement in Dicke states is highly resilient against 



external perturbations and measurements on individual 
qubits [10, 11]. Through projective measurements on 
some of the qubits in the system, we can obtain states 
of different entanglement classes. Thus, the Dicke states 
can serve as a versatile resource for the preparation of 
multipartite entangled states. 

In addition to the technique used in [7], other meth- 
ods for the generation of Dicke states with trapped-ion 
qubits have been proposed [12, 13]. The basic features 
of our approach follow those proposals. We consider a 
chain of N ion qubits each initialized in the state 1 1) and 
collectively cooled to the ground state of motion for a 
particular mode. By addressing a single ion in the chain 
with a laser pulse tuned to the m*'' higher-frequency mo- 
tional sideband (blue sideband) [14] we produce the state 
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where we have labeled the motional Fock state with a 
subscript M. From \ip), a laser pulse of appropriate dura- 
tion addressing all N qubits and tuned to the first lower- 
frequency motional sideband (red sideband), creates a 
state that is a close approximation of [Z)^'). 

We define the fidelity for an arbitrary final state, p, as 
F EE (D<r'|p|D^'). In the case with m = 1 (W states), 
the procedure outlined above can achieve arbitrarily high 
fidelity. For higher-order Dicke states attainable fidelity 
is reduced because the red sideband pulse does not trans- 
fer population completely from to [D^vJ"'). The previ- 
ous proposals showed that this imperfection can be miti- 
gated by postselection [12], or by a generalization of adi- 
abatic rapid passage [13, 15]. However, the simplified 
method presented here is sufficient to achieve fidelities 
limited by other experimental imperfection for a range of 
Dicke states. 

One experimental challenge in the above scheme is to 
individually address a single ion to produce the state [t/i). 
High confinement frequencies, and consequently small 
inter-ion spacing, are desirable to resolve the motional 
sideband spectrum. This makes individual spatial ad- 
dressing difficult. We avoid this problem by introducing 
to the A'^-qubit ion chain a single ancillary ion of a differ- 
ent atomic species. Simultaneous trapping of individual 
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FIG. 1: (a) Fidelity of l-D^') state generation (by use of the 
in-phase axial mode of motion) as a function of mass ratio, 
fi, for N = 2 . . . 6. In the case of A'' = 2 we consider both 
the symmetric case labeled S, with the ancillary ion in the 
center, and the asymmetric case labeled A, with the ancillary 
ion on the outside. The latter case is the one chosen in the 
experiment. For N = 3to N = Gwe assume that the ancillary 
ion occupies a central position in the ion chain, (b) Fidelity 
of and I/^n') state generation as a function of mass 

ratio for A'^ = 4 . . . 6. Here the ancilla ion occupies a central 
position for all cases. 



ions of two species has been used for sympathetic cooling 
[16] and indirect state detection [17]. Spectroscopic res- 
olution of the two atomic species ensures that any laser 
pulse applied to one species will leave the internal states 
of the other ion species unchanged. The state can be 
created by applying a global laser pulse tuned to the m*'' 
sideband (or equivalently m pulses sequentially tuned to 
the first sideband) of an accessible transition in the an- 
cilla. 

Two-species ion chains typically contain ions of un- 
equal mass, which leads to unequal motional mode cou- 
plings. To see how this affects state fidelity, we consider 
the generation of W states ([D^')), for which the fidelity 
can be perfect in the case of equal masses. For the i*^ 
ion in the chain, we denote the ground-state motional 
amplitude as Zi. The Lamb-Dicke parameter is defined 
by rji = kzi, where k is the laser beam fc-vector pro- 
jection along the z direction. In the Lamb-Dicke limit 
{rji <C 1), the strength of the ion's couphng to the red 
sideband pulse is fti = Jlo^ij where flo is the carrier 
(|i)|"^)Af — > |T)l"^)A/) Rabi rate. Beginning from the 
state lip) with to = 1, under the red sideband interaction. 



the A'^-ion state acts like a two-level system as described 
by the Morris-Shore transformation [18]. The motion un- 
dergoes Rabi oscillations between the two Fock states |1) 
and |0) at a frequency, fl' that obeys the relation 
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in which the terms in the superposition are the same as 
those in but the state amplitudes are proportional 

to the individual ion coupling strengths, Hi. The obtain- 
able W-state fidelity, F = |(Z?^'|'(/'')|^, can be written in 
terms of the coupling strengths as 



F = 




(5) 



To quantify the effect on fidelity of different mass ra- 
tios and different ion configurations, we first calculate the 
motional amplitudes [19]. Each configuration is speci- 
fied by the number of qubit ions, N , and the mass ratio 
fj, = ManciiiaZ-Mqubit • With ouc cxccption notcd below 
we assume that the position of the ancilla is at the cen- 
ter of the chain (for N even) or adjacent to the center 
of the chain (for N odd). These configurations can be 
prepared deterministically by adjusting trap parameters. 
We calculate equilibrium positions for the -I- 1 ions, 
then determine the amplitude of small oscillations about 
equilibrium. The in-phase, axial mode of motion exhibits 
the smallest deviations between individual ion motional 
amplitudes and is used for all fidelity calculations. These 
amplitudes give us the coupling strengths, fli, and F fol- 
lows from Eq. (5). The results are presented in Fig. 1(a), 
where we have assumed the Lamb-Dicke limit for simplic- 
ity, although fidelities for {rji} « 1 are similar. The case 
labeled 13*^' S is the only one that yields a theoretically 
perfect fidelity for any ^ 1 because that symmetric con- 
figuration, with the ancilla in the center of two qubits, 
gives equal motional amplitudes for the two outer ions. 
The cases with A^ = 3 ... 6, as well as that labeled D'^^ A, 
with the ancilla at the outside position, allow high fidelity 
over a wide range of mass ratios. 

In Fig. 1 (b) we present the result of a similar calcu- 
lation for the states and [D^") with A^ = 4. . .6. 
During the application of the red sideband pulse the ions 
state evolves in a complicated manner. To determine a 
practical upper limit on the fidelity, we numerically find 
the first maximum of the quantity F — (D^'|p(t)|Z3^') 
as a function of the red sideband pulse duration. Here 
p{t) is the reduced density matrix of the qubit system at 
pulse duration t after tracing over the motional degrees 
of freedom. We find fidelities for the case m = 2 as high 
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as 0.99 and for m = 3 above 0.96. The effect of unequal 
masses is found to be similar to the W state case with a 
mass ratio of 10 reducing the optimum fidelity at the 1 % 
to 2 % level. These numbers refer to the fidelities at the 
first maximum during the red sideband evolution, but, in 
the absence of decoherence, F could in general be made 
higher by evolving for a longer duration and reaching a 
later maxima. 

Note that the fidelities for /i = 1 in Fig. 1 are equal to 
the fidelities obtainable in the case where there is no an- 
cilla and is generated by spatially addressing a single 
ion with a strongly focussed laser beam. Here, there are 
only N ions in the chain, and they share equal motional 
amplitudes for the in-phase axial mode. Under the red 
sideband interaction beginning from j'i/;) the system can 
be transformed to a basis where it evolves as a ladder of 
m + 1 states spaced equally in energy [20] . Here the red 
sideband Hamiltonian is symmetric with respect to ion 
exchange, so all ions participate equally in the entangled 
state. 

We have demonstrated the basic features of the proto- 
col in an experiment with one ^^Al"*" ancillary ion and two 
25]y[g+ qubit ions. The ions are trapped in a linear RF 
Paul trap [21], with trap frequencies for a single ^^Mg+ 
ion {ujx,ujy,u;z} = 27r x {5.54,6.46,2.55} MHz where z 
denotes the axis of the ion chain. In ^'^Mg+ we use two 
Zeeman sub-levels of the ground-state hyperfine manifold 
as a qubit. We define |J.) = |^Si,F = 3, nip = —3) and 

It) = I^Si , F = 2, TTiF = -2). In 27A1+ the relevant qubit 
levels are defined as l^^o) = |^So,F — ^,mp — — |) and 
|3Pi) = |3Pi,F = |,TOF = -|). The ions are loaded 
into the trap via photo-ionization. 

The order of the ions is maintained as Mg-Mg-Al by 
monitoring the motional spectrum and adjusting DC 
voltages of the trap electrodes to regain the correct or- 
der when necessary. Specifically, we raise the DC endcap 
voltages and apply a radial bias field to configure a radi- 
ally oriented, linear ion chain with Al at one end. Then 
we apply differential endcap voltages to twist the radial 
chain to the desired orientation and finally relax the volt- 
ages back to the experimental parameters. In the order 
Mg-Mg-Al, the amplitudes of motion for the Mg+ ion in 
the in-phase mode are equal to within 1%, having a neg- 
ligible impact on fidelity, while the Mg+ ion spacing is 
small (3 /im) to facilitate equally strong interaction with 
the laser beams. 

Initially, the three axial motional modes of the ions 
are cooled close to the ground state by resolved side- 
band Raman cooling of Mg+ [22]. This process fills 
the largest part of our experimental duty cycle, about 
2 ms. We observe residual phonon numbers n < 0.1 
for all axial modes. Optical pumping ideally prepares 
the system in 1 11)1^50)10) j\f. The preparation of \D^^^) 
starts with a laser pulse that adds one phonon to the in- 
phase axial mode by driving the blue sideband transition 
r'S'o)|0)M |"^-Pi)|1)m in the A1+ ion, changing the ion 
state to I ii)|^-Pi)|l)Af- Then a laser pulse removes one 
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FIG. 2: (a) Example experimental histogram of photon counts 
from fluorescence detection of entangled ions (bars) and its 
fit (dot-dashed line). The stepped lines represent calibrated 
histograms P{n\0), P(n|l), and P{n\2), scaled respectively 
by the fit parameters {co,ci,C2} = {0.08,0.80,0.12}. The fit 
gives a parity of 11(0, 0) = -0.80(1). (b) Parity of qubit state 
after applying a rotation R{^,(f)) of varying phase (squares) 
and the parity after first applying the rotation i?(-|, §■), then 
applying a second rotation, i?(-|, 0), of varying phase (circles). 
The residual oscillation in the former case (squares) arises 
from small even-parity populations in p. The measurements 
give a an overall fldelity F — 0.77(2). 



phonon from the in-phase mode by driving the red side- 
band transition in the Mg"*" ions. This ideally produces 
the state ^(| {{) + | T}))I'A)|0)m [23]. 

In terms of the individual density matrix elements, the 
fidelity of the final Mg+ state, p, is 

^ = 2(''J-TaT + Pu-.n + Pihu + Pn.n)- (6) 

We measure the fidelity with the same technique used 
in previous experiments [6]. The odd-parity population, 
PiT iT + PtbTi' measured using resonance fiuorescence 
histograms from p. The off-diagonal term, Pj,t,Ti, +PTi>iT' 
is inferred from the parities of the states obtained by 
applying a qubit rotation R{^,(f)) of variable phase. 
Here we use the convention that R{0, (p) implements the 
transformation, ||) cos(|) ||) — ze~"^sin(|) [{) and 
It) -ie+*'^sin(f ) [}) -I- cos(f ) |t). In terms of the par- 
ity operator for two qubits, 

n = iii)aii + itt)(tti- (UT)aTi + iTi)(ni), (7) 

we define H(6', </>) = tr (i?t(6», cfypRid, 0)n) . We can then 
write px|,tx+Pn,iT = 5^(1, 0) -f n(f , f )]. A signature 
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of the state l-D*"^') is that it produces an even parity state 
independent of the 7r/2 analysis pulse phase. 

We detect the state of the Mg+ qubit by applying 
cr— polarized laser pulse of duration 200 /xs resonant with 
the IJ,) pP3,F — 4,m_F = —4) cycling transition and 
counting photons. We fit observed photon histograms to 
a weighted sum of the distributions P{n\0), P{n\l) and 
P{n\2), which correspond to the probability of observing 
n photons given 0, 1 or 2 ions in the state |J,). Assuming 
equal illumination and equal photon collection efficiency 
for the two ions, we can write the distributions above in 
terms of single-ion count distributions as 

P{n\0) = PBG{n) * Pt W * -Pt W («) 

P{n\l) = PBGin) * P^{n) * P^{n) (9) 
P{n\2) = PsGin) * Pi(n) * P^in), (10) 

where G{n) * = X^mXn C!{n — m)H{'m) is the dis- 

crete convolution of the distributions G and H. Here the 
convolved distributions refer to photon count probabili- 
ties from background scattering {BG) , as well as the two 
qubit states |J,) and |t). In the case of P]{n) we account 
for off-resonant repumping from |t) to |J,) by assuming 
an exponential decay rate, 7, such that if the ion be- 
gan in It), the state probabilities at time t are given by 
C| ~ e^''* and C| = 1 — C| [24]. The photon coimt rates 
and repump rate that determine the distributions result 
from a simultaneous fit to two reference histograms taken 
just before or after the experiment. 

In any particular measurement, we fit the unknown 
distribution 

2 

Pp{n) =Y,c,P{n\i) (11) 

4=0 

with ^ . Cj = 1 to the observed series of Np samples of 
photon counts, {ni, n2 . . . un^}- We use a maximum like- 
lihood method where we maximize the quantity 

L = \{Pp{n3) (12) 

by adjusting the parameters Cj. 

An example of the results of mcasiiring photon counts 
from the entangled state as well as a fit to the distribution 
are shown in Fig. 2(a). We measure odd-parity popula- 
tions of the initial entangled state of ci = —0.80(1). In 
Fig. 2(b) we display the measured parity of our state 
as a function of the phase of the analysis pulse (squares) 
as well as a sinusoidal fit. The mean value of 11(^,0) 
and n(|, |) is 0.74(2). These numbers together yield an 
overaU fidelity of F = 0.77(2). To further confirm the 



presence of entanglement, we applied the pulse i?(^,0) 
to p followed by a second analysis pulse i?(.|, 0) of vary- 
ing phase (circles), again measuring the parity of the final 
state. Here we observe quantum state interference in the 
sinusoidal parity oscillation with a period of tt and an 
amplitude of 0.70(3). 

Several experimental imperfections contribute to the 
overah infideHty oi 1 - F = 0.23(2). The largest er- 
ror stems from random failures to optically pump to the 
A1+ \^So,mF = — |) ground state at the 10 % level, 
which has been tested in separate experiments. The next 
largest error is due to beam pointing fluctuations. If the 
two Mg"*" ions are unequally illuminated by the Raman 
beams during the entangling pulse, two errors arise. The 
state amplitudes will accumulate a differential phase due 
to unequal Stark shifts and the state probabilities will 
be different due to unequal coupling strengths. We bal- 
ance the Stark shifts by suppressing the beating of fringes 
in a Ramsey experiment, in which each Raman beam is 
applied separately during the wait period. Even at the 
optimum beam position, beam pointing fluctuations can 
cause errors, which we estimate to be approximately 5 %. 
Additional errors are caused by imperfect ground-state 
cooling of the in-phase mode of motion (0.05 quanta re- 
maining, 3 % infidelity), Al"*" decoherence from laser and 
magnetic field noise (3 % infidelity), and A1+ sponta- 
neous emission (300 us upper state lifetime) during the 
|^'S'o)|0)m P-Pi)|1)m pulse {15 fis pulse duration, 1 % 
infidelity). These estimated errors add to 0.22 and agree 
with the experimental fidelity. 

In summary, we have explored numerically and exper- 
imentally a protocol for the creation of Dicke states in a 
trapped ion chain. The pulse sequence we use involves 
just two consecutive laser pulses, tuned to the resonance 
of a motional sideband and addressing all ions simultane- 
ously. Infidelities due to unequal ion masses and inequiv- 
alent positions in the chain can be small (<C 1%) for cre- 
ating W states, and the infidelity for creating Dicke states 
of two and three excitations can be sufficiently low to 
enable future interesting experiments. We have demon- 
strated the basic features of the process on an ion chain 
composed of two Mg+ ions coupled to a single A1+. If 
the technical errors observed in this demonstration are re- 
duced, scaling the experiment to a larger number of ions 
would require the same number of steps, which makes it 
an attractive method for enabling the study of multipar- 
tite entangled states. 
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ful comments on the manuscript and acknowledge the 
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Quantum Information program. This Letter is a contri- 
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